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We study a Ginzburg-Landau model of structural phase transition in two dimensions, in which a 
single order parameter is coupled to the tetragonal and dilational strains. Such elastic coupling terms 
in the free energy much affect the phase transition behavior particularly near the tricriticality. A 
characteristic feature is appearance of intermediate states, where the ordered and disordered regions 
coexist on mesoscopic scales in nearly steady states in a temperature window. The window width 
increases with increasing the strength of the dilational coupling. It arises from freezing of phase 
ordering in inhomogeneous strains. No impurity mechanism is involved. We present a simple theory 
of the intermediate states to produce phase diagrams consistent with simulation results. 



I. INTRODUCTION 

At structural phase transitions, anisotropically de- 
formed domains of the low-temperature ordered phase 
emerge in the high-temperature disordered phase, as the 
temperature is lowered [Ull]- In these processes, elastic 
strains are induced around the domains, radically influ- 
encing the phase transition behavior. The order parame- 
ter can be the concentration, the atomic configuration in 
a unit cell, the electric polarization, the magnetization, 
the orbital order, etc. Thus the effects are extremely 
varied and complex in real solids and have not yet been 
well understood. Among them, particularly remarkable 
are the so-called precursor phenomena taking place in a 
temperature window at first-order structural phase tran- 
sitions [3 H |5J H H [HI i]. They are caused by equihb- 
rium or metastable coexistence of disordered and ordered 
phases on mesoscopic scales and the volume fraction of 
the ordered regions increaes with lowering of the tem- 
perature, as unambiguously observed with various ex- 
perimental methods. Theoretically, phase ordering phe- 
nomena in solids have been studied using time-dependent 
Ginzburg-Landau or phase field models, in which an or- 
der parameter is coupled to the elastic field in a coarse- 
grained free energy functional [TJ [5] . Simulations of the 
model dynamics have been powerful in understanding the 
formation of mesoscopic domain structures. 

Kartha et al. (10] have ascribed the origin of the ob- 
served tweed patterns [4] to quenched disorder imposed 
by the compositional randomness. In their model, the 
elastic modulus C = (Cn — Ci2)/2 for the tetragonal 
deformations consists of a mean value proportional to 
T — Tq and a space-dependent random noise, where Tq is 
the nominal transition temperature. This randomness 
then produces quasi-static filuctuations of the strains, 
but it remains unclear how such a microscopic pertur- 
bation can lead to mesoscopic domain structures and 
lattice distortions above Tq. On the other hand, Seto 
et al. [6 proposed an intrinsic (impurity-free) pinnng 
mechanism stemming from anharmonic elasticity. One 
of the present authors demonstrated that third order 
anharmonic elasticity can freeze tetragonal domains in 



a disordered matrix in two-dimensional (2D) simulation 
[21 [TT]. Moreover, some authors stressed relevance of the 
anisotropic lattice deformations in producing multi-phase 
coexistence in perovskite manganites (possibly together 
with the compositional randoness) [T2J [T3]. For exam- 
ple, Ahn et al. [T3] have used a 2D phenomenological 
model with a two-component order parameter coupled 
to the tetragonal strain, but their model free energy con- 
tains only terms even with respect to the order parameter 
(improper coupling). As a result, in their model, the do- 
mains in phase ordering grow anisotropically up to the 
system size without pinning. 

In phase transitions subject to elsticity, there are var- 
ious intrinsic pinning mechanisms [5]. As such exam- 
ples, we mention binary alloys [TJ] and polymer gels 
[15|, where the elastic moduli depend on the composi- 
tion (elastic inhomogeneity) serving to freeze the domain 
growth. In gel- forming polymeric systems, furthermore, 
randomness in the crosslinkage constitutes quenched dis- 
order. Although the domain pinning itself can be induced 
by the elastic inhomogeneity only, the crosslink disorder 
produces quasi-static composition fiuctuations enhanced 
toward the volume-phase and sol-gel transitions [5] as de- 
tected by scattering experiments [TBI HZ] • We also men- 
tion hexagonal-to-orthorhombic transformations [21 [18] , 
where the interfaces between the ordered variants take 
preferred orientations hexagonal with respect to the prin- 
cipal lattice axes. This geometrical constraint gives rise 
to pinning of the domains. 

In this paper, we will further investigate the impurity- 
free pinning mechanism at structural phase transitions 
in the Ginzburg-Landau scheme. We will present ana- 
lytic results supported by extensive numerical calcula- 
tions. We will present only 2D results for the mathe- 
matical simplicity. We suppose solid films undergoing a 
square-to-rectangle transition in the plane. Such phase 
transitions would be realized under uniaxial compression 
[21 [19]. However, in real epitaxial films, analyzing elastic 
effects is difficult, because the displacement is fixed at 
the film-substrate boundary and the stress is free at the 
film-air interface [2U] . 

The organization of this paper is as follows. In Sec- 
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tion 2, we will set up the Ginzburg-Landau free energy 
and the dynamic equation in 2D. The elastic field will 
be expressed in terms of tp under the mechanical equi- 
librium condition. In Section 3, we will present analytic 
theory on two kinds of intermediate states at fixed vol- 
ume, in which ordered domains appear in a disordered 
matrix and take shapes of long stripes or lozenges. In 
Section 4, we will give simulation results in 2D under the 
fixed- volume condition, which will be compared with our 
theory. In Appendix B, we will derive the effective free 
energy in the strain-only theory. 



II. THEORETICAL BACKGROUND 

For ferroelastic transitions, one theoretical approach 
has been to set up a nonlinear elastic free energy contain- 
ing the strains up to sixth-order terms (the strain-only 
theory) [101 HH IIH [HllSl [H] , where some strain com- 
ponents constitute a multi-component order parameter 
and the corresponding elastic modulus decreases toward 
the transition. Another approach has been to intoduce 
a true order parameter different from the strains. If 
tp is coupled to the tetragonal or shear strains properly 
or in the bilinear form in the free energy, elastic soften- 
ing follows above the transition and anisotropic domains 
appear in phase ordering. For improper structural phase 
transitions [5], on the other hand, the square of the or- 
der parameter -0^ is coupled to the strains without elastic 
softening. Thus the second approach has been used for 
improper structural phase transitions [HI [25]. In this 
paper, we will adopt the second approach in 2D in the 
presence of both proper and improper elastic couplings, 
while one of the present authors presented a strain-only 
theory to describe intermediate states in 2D [TT] . 



A. Ginzburg-Landau model in 2D 

We assume that the free energy function F — F{ip, u} 
depends on a one-component order parameter ip and the 
displacement u. We define u as the elastic displacement 
from the atomic position in a reference one-phase state 
in the disordered matrix. No dislocation will be assumed 
in the coherent condition [2], where u is treated as a 
continuous variable. 

The free energy F consists of four parts. 



F 



dr 



/o + ?|VVP + /ci + /i„t 



(2.1) 



where the space integral is within the solid. The free 
energy density /o = /o('0) is the chemical part of ijj given 
in the Landau expansion form. 



/o 
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where t depends on the temperature T as 

T = Ao(r-ro). 



(2.2) 



(2.3) 



We retain the terms up to the sixth order in /q. The 
coefficients Aq, v, and C are positive constants, and Tq 
is the critical temperature in the absence of the elastic 
coupling (for m > 0). We call r the reduced tempera- 
ture. The other coefficients are assumed to be indepen- 
dent of T. For the mathematical simplicity, we assume 
the isotropic elasticity with homogeneous elastic moduli. 
Then the elastic energy density is of the harmonic form. 



/el = ye 



(2.4) 



where K is the bulk modulus and is the shear modulus. 
The ei, e2, and 63 are the dilational, tetragonal, and 
shear strains, respectively. In terms of the displacement 
vector u , they are defined by 



ei = "^xUx + "^yUy^ 
63 = VyUx + VxUy, 



(2.5) 



where = d/dx and Vy = d/dy. The third term J-„ 
consists of the following two interaction terms. 



fn 



-ge2ip — aeiip 



(2.6) 



Thus tp is properly coupled to the tetragonal strain 62 and 
improperly coupled to the dilational strain ei. Without 
loss of mathematical generality we may assume 



5 > 0, a > 0. 



(2.7) 



If both are negative, we change u to —it. For g < Q and 
a > we change ■0 to —tp, while for 5 > and a < the 
signs of u and "0 are both reversed. 

In the absence of the elastic coupling [a — g = Q) 
[3 [55], it is well-known that the transition is second-order 
for u > and is first-order for u < with the tricritical 
point at r = S = 0. In our problem, the phase transition 
behavior is much altered by /int ■ The coupling to 62 {g ^ 
0) gives rise to anisotropic domains (even in the isotropic 
elasticity) and that to ei (a ^ 0) makes the effective 
reduced temperature(= r — 2aei) inhomogeneous in the 
intermediate states [5]. 

We introduce the elastic stress tensor by 



'yy 



'xy 



= Kei + ^62 — aijj — gip, 
= Kei — ^62 — + gipj 



'yx 



(163. 



(2.8) 
(2.9) 
(2.10) 



If we change u by a small amount Su, the free energy F 
changes by SF — J drJ2ij'^ij9Suj/6xj. We then have 
SF/Sui = — Vjcr.y , where ip is fixed in the functional 
derivative with respect to m. In this paper, u is deter- 
mined from the mechanical equilibrium condition, 



SF 
6ui 



0. 



(2.11) 
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Then u becomes a functional of ■0 under given boundary 
conditions [HE]- 

In dynamics, the elastic field is assumed to instanta- 
neously satisfy the mechanical condition (2.11), while the 
order parameter ip obeys the relaxation equation. 



Ik 



-An 



5F 



(2.12) 



The kinetic coefficient Aq is assumed to be a constant. 
From Eqs.(2,2) and (2.6) /t in Eq.(2.12) is written as 

ji = Ttl) + ui})'^ + vi)^ - ge2 -2aeii}}. (2.13) 

If we integrate Eq.(2.12), equilibrium or metastable 
states can be reached at long times, where it holds the 
extremum condition, 



5F_ 



= fi- cv V = 0. 



(2.14) 



The gradient term is important in the interface regions, 
while we may neglect it and set fi = outside them. 



B. Elimination of the elastic field 

We may apply an average homogeneous strain as 
(VjMi) — Aij, where Aij is a homogeneous tensor repre- 
senting affine deformation. Hereafter (•••) = / dr{- ■ ■)/V 
represents the space average in the solid, V being the 
solid volume. The displacement may be divided into the 
average and the deviation as 



Ui 



7 . ^ij^j 



Sui. 



(2.15) 



The simplest boundary condition is to impose the peri- 
odicity of the deviation Su. Note that the space aver- 
age of J2ij crij\7iSuj vanishes in this case. Since 2/ei — 
o'ij'^iUj + aeiip'^ + ge2'0, the total free energy F in 
Eq.(2.1) becomes 



F ^ dr 
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Fr,. 



(2.16) 



Under the mechanical equilibrium condition (2.11), the 
elastic part F^ is written as 



F, 



dr 



(2.17) 



In this paper, we impose the periodic boundary con- 
dition on Su and tp in the region < x < Y^^"^ and < 
y < V^l"^ in 2D to investigate the domain morphology 
not affected by the boundaries. It is then convenient to 
use the Fourier transformation, 6uj(r) ~ l^k^jk^^^^ ^ 
where j — x and y in 2D. The mechanical equilibrium 
condition (2.11) may be rewritten in terms of u^f^ as 



where e^f^ and are the Fourier components of ei and 
■0, respectively, and is that of 



(2.19) 



Here we make (1^9) = 0. The Fourier components of the 
strains ei, 62, and 63 are calculated for fc ^ as 



Le^f^ = a(pfg + gtpf^ cos 29, 



^^2k ~ ^'^k 2^ + gi'ff 
K 



1 



K 



sin-^ 261 



Le. 



3fc 



sin 29 



^Vk 



A* 



gVfc cos 29 



where L is the longitudinal modulus defined by 

L = K + ^1. 
We set cos9 = kx/k and sin^? = ky/k, so 

cos 26*= {kl-kl)/k^, sm29 = 2kxky/k^, 



(2.20) 
(2.21) 

(2.22) 

(2.23) 
(2.24) 



in Eqs.(2.20)-(2.22). Here it is convenient to introduce a 
variable x('") defined by 



V^X = (V?-V^)0 



(2.25) 



with (x) = 0. The Fourier component of x satisfies Xk ~ 
i/jfgCos29 for fe 7^ 0. From Eq.(2.20) ei is expressed as 

ei = {ei) + {a^ + gx)/L. (2.26) 

In terms of x and (p, F^ in Eq.(2.17) is expressed as 



F,. 



dr 
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(e2)0-^(ei)02 
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(2.27) 



where Sip = "ip ~ {%p) . Thus F^ consists of contributions 
of the first, second, third, and fourth orders with respect 
to ip, where the coefficients depend on the space averages 
{ip) and (0^). 



C. Dimensionless representation 

In this paper, we will show snapshots of ^p and the 
strains obtained in our simulations. They will be nor- 
malized as Tp/ipo, ei/eo, and 62/69, where ipo and cq are 
typical amplitudes given by 

^y\u\K K y/\u\K^ 

in terms of the coefficients g, u, and K in the free energy 
F. We will set fi — K and include the case of u < 0. In 
terms of the coefficient C of the gradient free energy, the 
typical spatial scale i is given by 



(2.29) 



FIG. 1: Normalized order parameter 'ip/tpo in a stripe inter- 
mediate state for (t/tq, q/qq, 'w/^o) = (0.6,1.4,1). In the 
domains ^ — ips, —ips, or 0. 



FIG. 3: Normalized order parameter i/'/i/'o in a lozenge in- 
termediate state for {t /tq, a/ao,v/vo) = (0.6,1,1). In the 
domains tp = ^pe, —"ipt, or 0. 
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FIG. 2: Disordered square regions are enclosed by ordered 
stripe domains uniaxially stetched along the y or x axis. 



III. PHASE TRANSITIONS AT FIXED 
VOLUME AND SHAPE 

The phase transitions in sohds can crucially depend on 
the boundary condition |21 HZ] . In this paper, we limit 
ourselves to the simplest case of fixed volume (area in 
2D) and shape. This condition holds in the lateral direc- 
tions if a solid film is fixed to a substrate. For 3D solids, 
the phase transition is usually observed under the stress- 
free boundary condition and claming of the boundaries 
is needed to realize the fixed-volume condition. 

In this paper, we thus assume Aij = (i, j = x, y) and 



(ei) = (62) = (eg) - 0, 



(3.1) 



Furthermore, in our 2D simulation, there was no macro- 
scopic order or 



(^) = 0. 



(3.2) 



These expressions give the typical free energy density 
C^/'2/^2 ^ |^|^4 ^ g^^^^ ^ ^ £,4/1^1^2. We scale 

the coefficients t, a, and v by tq, ao, and vq, respectively, 
defined by 



9 

-0 = ^, 



(2.30) 



Our simulation results will be parametrized by t/tq, 
a/ao, and v/vq. If a ~ ao and -0 ~ "00; the ratio of 
the typical dilational strain (~ aip'^/K) to the typical 
tetragonal strain (~ gip/fJ-) is of order fi/K(= 1 in our 
simulation). Near the tricritical point, u tends to zero 
and ag and Vq become small so that the scaled parame- 
ters a/ao and v/vq are amplified. 



Then (aoox) = (cryy) = -a{'>p'^) and {a^y) = 0. Note 
that (ip) becomes nonvanishing under applied tetragonal 
strain (62) 7^ 0, since (62) plays a role of the ordering 
field in Eq.(2.27). In phase ordering, the condition (3.2) 
means that domains with positive ip and those with neg- 
ative ip equally appear and the space average of ■0 over 
many domains vanishes. Under Eqs.(3.1) and (3.2) the 
elastic contribution in Eq.(2.17) becomes 



F,= dr 



w 

2^^ 



9^ 2 



ag 



2L 



(3.3) 



where and x are defined by Eqs. (2.19) and (2.25), 
respectively. 

In our previous paper |27j . the model with a > and 
f; = has been studied at fixed volume. There, in a tem- 
perature window, the free energy is lower in two-phase 
states than in one-phase states, where the domains attain 



FIG. 4: Strains ei, 62, and F = ei — aip/K in units of eo in Eq.(2.28) in the same lozenge state as in Fig. 3. These stains are 
flat and F is smafl within the domains. As a result, ei = aip/K, 62 = g'<p/fJ- in accord with Eq.(3.35). 



macroscopic sizes, however. For a = and g > 0, on the 
other hand, lamellar or twin ordered states appear at low 
temperatures, where the interface normals make angles 
of ±7r/4 wth respect to the x axis and x vanishes in 2D 

[HI HI 112,^^^ ]. 



A. Stripe intermediate states 

Figure 1 is an example of intermediate states consist- 
ing of ordered stripes parallel to the x or y axis. As illus- 
trated in Fig. 2, the stripes enclose disordered rectangu- 
lar regions with ip = isotropically dilated at ei = e°. A 
rotation of 7r/2 of a stripe pattern yields another possi- 
ble stripe pattern. As a result, the areal fractions of the 
stripes along the x axis and those along the y axis are 
the same under no applied tetragonal strain. We assume 
that the width of the stripes is much longer than the in- 
terface thickness. Within the stripes along the y axis, we 
have 



62 



(3.4) 



On the other hand, within the stripes along the x axis, ip 
and 62 are reversed in sign, but ei is unchanged, so that 



62 



(3.5) 



Thus we have (ip) — and (62) — (63) = on the average. 
From the dilation-free condition (ei) = 6^-1- 4>es = in 
Eq.(3.1), the areal fraction of the ordered regions (j) is 
expressed as 

= -e;/e,. (3.6) 

To calculate the excess dilation Cg, we consider the in- 
terface region of a stripe along the y axis, where all the 
quantities depend only on x. The mechanical equilibrium 
condition yields a^.^. ^constant along the x axis. Here 
\7xUx changes along the x axis, while VyUj, remains at 
65/2. From Eq.(2.8) or from Eq.(2.26) VxUx across the 
interface is expressed as 



1 



61 



ip{x) 



L 



ip{x), 



(3.7) 



where L is defined in Eq.(2.23). The above quantity is 
equal to ei{x) — e\ and 62(0;). The result (3.7) also follows 
from Eqs.(2.20) and (2.21), where cos 20 = 1 and x = V' 
since ip change along the x axis. The order parameter 
i^{x) changes from to tl)g with increasing x, leading to 



(3.8) 



In the extremum condition (2.14) for %l^{x), we set ei — e\ 
and 62 equal to the right hand side of Eq.(3.7) to obtain 

= " ^""i"^' " + 2«V)(a^' + 5^)/i, (3.9) 

where /g = dfo/dip. This equation is integrated to give 
the standard form of the interface equation [5] , 



'2\dx ~ 



(3.10) 



where /st(V') is the effective free energy density near the 
interface. Here the derivative dfst (V') / dip is equal to the 
right hand side of Eq.(3.9) and /st(0) = is required, so 
we obtain 



fst = /o-ae?^'-(g + a^)W2i 



4" 6' L 
In the second line, we define the coefficients 



T 

u 



2a^/L, 
u-2a^/L- 



/3. 



L- 



(3.11) 



(3.12) 
(3.13) 
(3.14) 



The integrand of in Eq.(2.27) yields the above fst if 
we replace (ei) by ej, S^p by ip, x by V'l and ip by V'^, 
with (62) = 0. 

For the existence of a planar interface /st {ip) should be 
minimized at the two points "0 = and ips, so we need 
to require fst = dfst /dip = 'At ijj — ipg. Hence ips is the 
solution of the following cubic equation. 



d_fst 
dip ip^ 



2v 



.3 = 0. (3.15) 



1p = tlls 
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where r vanishes and ips turns out to be independent of 
T or the temperature. In Appendix A, we will explicitly 
solve the above cubic equation. From /stlV's) = the 
effective reduced temperature r' is related to ips by 



2a 
~L 



(3.16) 



where the quadratic term in the first line is eliminated 
in the second line with the aid of Eq.(3.15). Thus t' is 
also a constant independent of t. With Eqs.(3.13) and 
(3.14), we may rewrite fst as 



, , , ^2 " ,2 W 



(3.17) 



Since fstii') > should hold ai ip — —ips, the inequality 
vips/S + u/4: > follows, which then gives ips > from 
Eq.(3.15). The positivity r' > also holds. 

From Eq. (3.12) we have (r - g'^/L - T')/2a. 

Use of Eqs. (3.8) and (3.16) gives the areal fraction (p in 
Eq.(3.6) in the form. 



= (Ts - t)/Ts, 



(3.18) 



Here = t' + /L and Tgy^ = 2aes are rewritten as 



3^- 



q a 
2 

Tsw = -aipsig + aips) 



(3.19) 
(3.20) 



Since (/) — > as r — > Tg, Tg is the transition-point value of 
T. From < ^ < 1, T is in the window range. 



(3.21) 



to realize the stripe domain morphology. For r <Ts — 
Tsw the system changes into an ordered state (where the 
order parameter and the free energy will be calculated as 
Eqs. (3.47) and (3.48), respectively). 

In the stripe intermediate states, the free energy den- 
sity in the disordered regions is fei — K{e\)'^ /2, while 
that in the stripes is /o -I- /ei + /int — Kei{ei/2 + Cs) 
from /st(V's) = or from Eq.(3.16), where Cg is the ex- 
cess dilation given in Eq.(3.8). Neglecting the surface 
contribution, we write the free energy density as 



F 
V 



K 



K 

'8a 



2(Ts 



2(pes) 
-r)\ 



(3.22) 



where the first line holds for general and the second 
line follows from Eq.(3.6) at fixed volume. Thus F is 
negative in the stripe intermediate phase. 




FIG. 5: Illustration of the lozenge geometry. Symmetric pat- 
tern consisting of lozenges and squares (left) and asymmetric 
pattern consisting of parallelograms and rectangles (right). 
The angles of the interface normals with respect to the hori- 
zontal axis are written in terms of the angle 5. The numbers 
1, —1, and denote the value of ip{x, y)/ipi. 
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FIG. 6: Normalized order parameter ip/ipo in a twin state 
consisting of two ordered variants separated by interfaces for 
{T/To,a/ao,v/vo,) = (—1,1,1), where ipo, tq, ao, and vo are 
given in Eqs. (2.28) and (2.30) and the spatial scale is I in 
Eq.(2.29). 



B. Lozenge intermediate states 

We show another kind of intermediate states. In Fig, 3, 
the pattern consists of ordered lozenges (rhombus) with 
ip — zLipi and disordered squares with -0 = 0. As in Fig. 
4, the strains ei and 62 change mostly in the interface 
regions and are nearly fiat (homogeneous) within each 
domain. The geometry is illustrated in Fig. 5. The or- 
dered domains are characterized by the lozenge angle 5 
in the range < S < tt/A. More generally, as shown on 
the right, the pattern may consist of ordered parallelo- 
grams (quadrilaterals with both pairs of opposite sides 
being parallel) and disordered rectangles. The stepwise 
behavior of the strains in each domain in Fig. 4 is a char- 
acteristic feature of the lozenge geometry in Fig. 5. We 
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FIG. 7: Reduced transition temperature Ttr (maximum of Ts, 
T«, and roc) in units of To as functions of a/ao. The numbers 

to each curve represent (u/\u\,v /v^). The transition is from 
disordered to lozenge states on the sohd lines (L). At the 
upward arrow 5 7r/4, and on its left (O) the transition is 
to twin ordered states. At the downward arrows S —^ 0, and 
on their right (S) the transition is to stripe states. 



may derive it starting with the assumption of the step- 
wise behavior of tp. We first examine the behavior of x in 
Eq.(2.25) in the vicinity of an interface between ordered 
and disordered domains. Except near the corners of the 
lozenges, we may set 



^(a;,t/) =±V^*(C), 



(3.23) 



where + (or — ) is for the ordered phase with tjj = tpi 
(or Tp = —^i) and C is the coordinate along the in- 
terface normal direction n ~ (n^^ny). The function 
\E'(C) tends to unity in the interior of all the ordered do- 
mains, while it vanishes in the disordered regions. From 
Eq.(2.25) we obtain (fx/dC"^ = ±cos{2eo)i^ecP'if /dC 
where cos(2^o) = i^x~''^y "^^'i^^^ being the angle between 
n and the x axis. From Fig. 5 we can see that = 
(or 6*0 = 7r/2 ± (5) for the ordered phase with ip ^ tpi 
(or for that with tp = ~tpe) so that cos(20o) = ±cos(25). 
Thus we find 



d^ 



cos{26)tpi 



' dC '■ 



(3.24) 



where the symbol ± has disappeared. This is satisfied if 
X — cos{25)ipe,^ is a constant or if 



X{x,y) = [^{x,y) - 4>]'il)tcos{25), 



(3.25) 



where cf) = (^') is the areal fraction of the ordered regions 
and (x) = is satisfied (since \E'(a;, y) is equal to 1 in the 
ordered regions and to in the disordered regions). We 
numerically calculated x(x, y) from the inverse Fourier 
transformation of Xk — i^x ~ ^y)^k/^^ *° confirm the 



above relation for a number of the lozenge states, where 
the agreement becomes better with decreasing the inter- 
face width (or decreasing the coefficient C in Eq.(2.1)). 
Similarly, the strain e2{x,y) is proportional to 'ip{x,y) 
and is of the form, 

e2(x, y)^[g- KT, cos(25)] ^(a;, y)/^. (3.26) 

where we define 

Fe = [g cos(2<5) - a/^V^/i^] V'^/^- (3.27) 

Here the inverse Fouri(^r transformation of oup^^ cos(20) in 
Eq.(2.,21) has been approximated by acos(25)tpitp{x,y). 
In the thin interface limit, i,b{x,y) /tpf {= 0, ±1) and 
^'(x,j/)(= 0,1) are step functions and Eqs.(3.25) and 
(3.26) hold exactly as solutions of Eqs.(2.21) and (2.25). 
We rewrite the elastic free energy in Eq.(3.3) as 



dr 



2/x 



2/x 



(3.28) 



where we define the coefficient Pk and introduce the 
quantity F by 



Pk = 2a''/K, 
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-X- 



if. 



(3.29) 
(3.30) 



Then ei F -|- aip/ K in terms of F. From Eq.(3.24) 
we find F = (1 — (j))Tc in the ordered regions and F = 
—4>Vf. in the disordered regions with Fj, being defined 
by Eq.(3.27), so that (F) = and (F^) = 0(1 - (j))Tl. 
In addition ((^^) = tpj(t){l — (p)- Neglecting the surface 
tension contribution, we may calculate the average free 
energy density as 



(1-0). (3.31) 



For simplicity, we introduce 



(3.32) 



which is of the same form as /o in Eq.(2.2) with r being 
shifted to r — g^Z/x. 

The average free energy in Eq.(3.31) depends on S, ijj£, 
and (p. In our simulation to follow, we shall see that both 
the symmetric and asymmetric patterns in Fig. 5 both 
appear and the edges of the ordered regions are flattened 
with decreasing r (Figs. 11 and 14). This means that 
S and (p can be varied independently, though they are 
related as (p = sin(2(5)/(l + sin(2(5)) for the symmetric 
pattern (left of Fig. 5). Thus we minimize F with respect 
to these three quantities. First, from Eq.(3.27) we notice 
that Fc can vanish for 



afi-ipe/gK < 1, 
if the angle S is chosen as 



5 = - cos ^{a/j.'ipi/gK). 



(3.33) 



(3.34) 
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If we set r = or X = {cttJ-/gK)ip, we obtain 



(3.35) 



where the former follows from Eq.(2.26) and the latter 
from Eq.(2.21) or from Eq.(3.36). The snapshots in Fig. 
4 are consistent with the above relations except in the 
interface regions. 

Under Eqs.(3.34) and (3.35) we furthermore minimize 
F/V with respect to tp^ and (p. If Tc = 0, the coefficient 
of the quadratic term (cx 0^) is positive in Eq.(3.31) and 
F/V can take a minimum as a function of (f> in the range 
<(/!>< 1 for sufficiently small f{ipi). From Eq.(3.31) 
the minimum conditions read 

1 dF 

-n7^ = /'(^£)-/3ifV'l(l-</') = 0, (3.36) 
(j)V dipt 

'\ f)j^ 1 

/(V'^) - - 20) = 0, (3.37) 

where /' = df /dip. The first equation (3.36) is equiva- 
lent to the equilibrium condition (2.14) outside the inter- 
face regions if use is made of Eq.(3.35). In our previous 
work on the compressible Ising model [57], similar mmin- 
imization of F yielded two-phase coexistence in a tem- 
perature window (leading to Eqs.(3.49)-(3.52) in the next 
subsection). In the present case, we may eliminate (p from 
Eqs.(3.36) and (3.37) to obtain i}jf - 2/ - fiK^^/2 = 
at ip = ipe. This is solved to give 



ipj = 3{f3K - u)/4v, 



(3.38) 



from which we need to require Pk > u. Remarkably, 
ip£ becomes independent of r as well as ips determined 
by Eq.(3.15). From Eq.(3.36) the areal fraction (p is ex- 
pressed as 



= 1 - f'{iJji)/Mi = in - T)/n^ 

in the same form as in Eq.(3.18). Here, 



(3.39) 



^- + li't=^- + ^{PK-ur, (3.40) 



f^Ki'j = -r/^KiPx - u). 



(3.41) 



The Ti is the transition reduced temperature and Tg^ is 
the width of the window. The lozenge intermediate phase 
is realized in the window, 



T£ — Tl^ < T < Tg. 



(3.42) 



As in Eq.(3.22) the average free energy density is calcu- 
lated as 



F 
V 



K 



{n-rf 



(3.43) 



Now we return to Eq. (3.33) and the inequality (3^ > u 
implied by Eq.(3.38). They impose the range of — 
2c? /K allowed for the lozenge phase. 



u<P3k < 




(3.44) 



Here Eq.(3.33) gives the upper bound, at which 8 tends 
to and the lozenge patterns continuously change into 
those of stripes. Therefore, the lozenge-stripe boundary 
is determined by ipi — gK / ^a, at which ■0s = ipe, leading 
to Ts — from Eqs.(3.19) and (3.40). The continuity of 
F/V also holds from Eqs.(3.22) and (3.43). 



C. One-dimensionally ordered states 

With decreasing t, our system eventually consists of 
regions of two ordered variants with ip = ±ipc at fixed 
volume and shape in steady states. In Fig. 6, we 
show such a lamellar or twin ordered state, where the 
two ordered variants are separated by antiphase bound- 
aries or twin walls and there is no bulk disordered region 
[HI [211 1221 [m m]. Here we assume that all the quanti- 
ties depend only on x' — {x — y)/^/2 as in Fig. 6. Then 
we obtain x = 0, leading to ei — aip/L, and 62 — gip/ 
The resultant free energy at fixed volume is expressed as 



F 



dr 



h 2 



(3.45) 



where ip' — dip/dx' and P3l is defined by Eq.(3.13). The 
above form coincides with the free energy of the com- 
pressible Ising model at fixed volume (with no ordering 
field conjugate to ip) ^7] , on the basis of which we discuss 
the two cases below. 

If u = u — /3i in Eq.(3.14) is positive, a twin ordered 
state is realized for 



(3.46) 



where we may set = and then ipc is the solution of 
f'{ipc) = with / being defined by Eq.(3.32). It is solved 
to give 



Av{T-g'^/li)Y/'^/2v~u/2v. (3.47) 



The elasticity effect is only to shift t to t — g^ / p.. The 
average free energy density takes a negative value. 



F/V ^-i,l{Zu + ^viPl)/\2, 



(3.48) 



where the surface tension contribution is neglected. 

For u — u — P3l < 0, on the other hand, the elastic term 
proportional to /3l in Eq.(3.45) can give rise to macro- 
scopic coexistence of disordered and ordered regions in 
the following window. 



(3.49) 



where ip — ±(3|u|/4ti)^/^ in the ordered regions. The 
transition-point value and the window width in this case 
are given by 

, (7^ 3u^ 3 , ^ n9 /„ 

roc - - + T7^ = - + T7r(^^-") ' ^^-^^^ 
/i low /i Ibv 

Ton, = ^\U\^ ^pLih-u). (3.51) 



9 



1.20 



1.02 



0.84 



0.66 



U > 0, v/l\) =0.1 
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FIG. 8: Theoretical phase diagrams in the plane of t/tq and q/qq for {u/\u\,v/vo) = (1,0.1) (left), (1,1) (middle), and 
(— 1, 1) (right). For tt > 0, we have the disordered (D), twin ordered (O), stripe intermediate (S), and lozenge intermediate (L) 
phases. For u < (right), there appears the lamellar intermediate phase (LI), where the disordered and ordered regions coexist 
macroscopically '2T. The points Ai, Bi, B2, Ci, C2, and C3 represent intermediate states with patterns shown in Figs. 9-13. 
At the points A2 and B3 twin patterns were realized numerically. 



The areal fraction of the ordered regions is cf) — (tq^ — 
t) /tqw ■ The average free energy without the surface ten- 
sion contribution becomes 



F 
V 



L 



(3.52) 



Below the lower bound of the window in Eq.(3.49) the 
twin ordered phase with ip — dzip^ is realized. In our 
simulation, however, we did not realize the macroscopic 
coexistence predicted above. 



D. Theoretical phase behavior 

As T is lowered from a value in a disordered state below 
a reduced transition temperature Ttr, phase ordering oc- 
curs into an intermediate or ordered state. Depending on 
which is largest, rtr is given by either of Tg in Eq.(3.19), T( 
in Eq.(3.40) under Eq.(3.44), tqc in Eq.(3.46) for u > 0, 
or t'^^ in Eq.(3.50) for m < 0. In Fig. 7, we show it as a 
function of a/ao for four sets of {u/\u\,v/vo) at 11 ~ K. 
(i) For M > 0, the transition is at t = j to the twin 
phase for a/ag < %/2, to the lozenge phase for the inter- 
mediate range. 



V2 < a/ao < Vl + A, 



(3.53) 



where A = {I + 32K'^v/3fj.'^vo)'^^'^ from Eq.(3.44). For 
larger a/ao the stripe phase is realized, (ii) For u < 0, 
the transition is to the lozenge phase for a/ao < + A 
and to the stripe phase for a/ao > Vl + A. 

We may calculate a phase diagram in the plane of 
t/to and a/ao using the scalings in Eqs.(2.28)-(2.30) if 
we give n/K, v/vo, and the sign of u. Figure 8 dis- 
plays such examples at {u/\u\,v/vo) — (1,0.1), (1,1), 
and (—1,1) at fi/K = 1. For u > 0, the plane is di- 
vided into the disordered (D), twin ordered (O), stripe 




/ 1-2, . 
a/ao 



FIG. 9: Theoretical order parameter tpM scaled by tpo (up- 
per plate), which is tps, i^e, or tpc in the stripe, lozenge, or 
twin phase, respectively. Theoretical areal fraction <f) of the 
ordered regions (lower panel) changing from to 1 as r is 
decreased. Here "u > and v/vo = 1 and the corresponding 
phase diagram is the middle plate of Fig. 8. 
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FIG. 10: Relaxation of the total perimeter length after 
quenching from a disordered to stripe state at Ci (lower curve) 
and Cs (upper curve) in Fig. 8. Evolution of tp is displayed. 
Times are in units of to in Eq.(4.1). 



FIG. 11: Relaxation from a disordered to lozenge state at 
Ai (lower curve) and Bi (upper curve) in Fig. 8. The final 
pattern is symmetric at Ai and asymmetric at Bi (see Fig. 
5). 



(S), and lozenge (O) regions, depending on which mor- 
phology has the lowest free energy (the interface free 
energy being neglected). The lozenge region can ap- 
pear for a/ao > from Eq.(3.44) and is narrow for 
small v/vq (or far from the tricriticality) . The boundary 
between the lozenge and twin phases is determined by 
T = Tg — Tiw (or (j) = 1)). In fact, at this temperature, 
Eq.(3.47) gives Ve = ^{Pk - u)/Av = V? and F/V in 
Eq.(3.48) becomes —Kt^/So^ , demonstrating the conti- 
nuity of ip and F/V at the boundary from Eqs.(3.38) and 
(3.52). For u < 0, the lamellar region is further divided 
into the ordered region (O) and the intermediate region 
of macroscopic two-phase coexistence (LI) described by 
Eqs.(3.49)-(3.52). The latter phase (LI) was not realized 
in our simulation, however. 

Furthermore, in Fig. 9, we plot the equilibrium or- 
der parameter ■f/'Af in the upper plate and the equilib- 
rium areal fraction (j) in the lower plate for u > and 
v/vq = 1. The former is given by -08, ^i, or tpc in the 
stripe, lozenge, or twin phase, respectively. It is contin- 
uous at the disorder-twin and lozenge-twin boundaries 
and discontinuous at the other phase boundaries. The 
(j) is between and 1 in the intermediate phases. It is 



continuous at the disorder-stripe, disorder-lozenge, and 
lozenge-twin phase boundaries and discontinuous at the 
other phase boundaries. We can see that ipM and cj) are 
both discontinuous at the twin-stripe boundary. 



IV. NUMERICAL RESULTS 

We numerically integrated the dynamic equation (2.12) 
in 2D under the periodic boundary condition on a 256 x 
256 lattice. We measure space and time in units of £ in 
Eq.(2.29) and 



to = f/\oC - K/Xog^ 



(4.1) 



where Aq is the kinetic coefficient. When we started 
with a disordered phase, the initial value of ■0 was a 
random number in the range {tpl < O.Olipo at each lat- 
tice point. No random source term was added to the 
dynamic equation. We calculated the strains ei and 62 
using their Fourier components in Eqs.(2.20) and (2.21) 
under the fixed volume condition. We set (ei) = (62) = 
0. We will show the resultant patterns at the points 
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FIG. 12: Numerical arcal fraction of the ordered regions (j) as a function of (r — rs)/Tsw (left) and (r — Ti)/Te^. (right) in the 
stripe and lozenge phases, respectively. The inirribers are the values of (a/an, v/v{)). See Fig. 8 for the positions of Ci, C2, C?,, 
Ai, A2, Bi, B2, and B3 in the phcise diagrams. The arrows in the left panel indicate how the data were obtained in the range 

5^//tt <T <Ts. 



marked by x in the phase diagrams in Fig. 8 (ex- 
cept A2 and B3). The parameters (t/tq, a/ao, w/wo) 
are given by Ai: (0.6,0.76,0.1), A2: (0,0.76,0.1), Bi: 
(0.6, 1, 1), B2: (0.4, 1, 1), B3: (-0.2, 1, 1), Ci: (0.6, 1.4, 1), 
C2: (-0.4, 1.4, 1), and C3: (-1.4, 1.4, 1). 

A. Transient behavior 

In transient states, the ordered domains with positive 
(negative) tp tend to be elongated along the horizontal x 
(vertical y) axis after their formation. They do not pen- 
etrate into the others belonging to the different variant 
(having the opposite ij)) on their encounters. Further- 
more, thickening of the stripe shape is prohibited by the 
fixed-volume condition (ei) = 0. As a result, our system 
tended to a nearly pinned state at long times in each 
simulation run. 

In Fig. 10, we display the time evolution of the to- 
tal perimeter length in the relaxation from disordered 
to stripe intermediate states at the points Ci (where 
t/to = 0.6) and C3(where t/tq = —1.4) in the mid- 
dle panel of Fig. 8. The initial stage of domain forma- 
tion is in the range t/U) <10 and the coarsening stops 

for t/tfy >10'^. At the higher temperature (Ci), the or- 
dered domains consist of long stripes with a smaller areal 
fraction. At the lower temperature (C3), the domains 
are distorted, consisting of both large and small ordered 
stripes. In Fig. 11, we display the relaxation from dis- 
ordered to lozenge intermediate states at the points Ai 
and Bi in the left and middle panels of Fig. 8. Here the 
domain pinning occurs at t/fo = 2 — 4 x 10^. The pat- 
terns are nearly symmetric at Ai but largely asymmetric 
at Bi. In these two cases, the lozenge angle 5 and the 



areal fraction (j) are nearly the same (see Fig. 12), while 
the domain sizes are considerably different. 



B. Steady intermediate states 

In our system, the disordered phase [ip = 0) is linearly 
stable for t > g'^/fi (as can be known from Eq.(3.3)). 
Therefore, we lowered t below g'^ / ^ to produce ordered 
domains. Notice that the intermediate phase can be sta- 
ble or have a negative free energy for r above g'^/fi and 
below the transition-point value, as can be seen in Figs. 
8 and 9. To realize such intermediate states, we increased 
T starting with twin or intermediate states. Then </> de- 
creased on approaching the transition point from below. 

In Fig. 12, we show <j) versus (t — Tg) /tsw in the steady 
stripe states (left) and <j) versus (r — re) /r^w in the steady 
lozenge states (right). Our theory predicts the linear de- 
pendence in Eqs.(3.18) and (3.39). The points Ai, A2, 
Bi, B2, B3, Ci, C2, and C3 correspond to those in the 
phase diagrams in Fig. 8. Slightly below the transition, 
these numerical data excellently agree with our theory. 
In addition, in the left panel, we show the hysteresis be- 
havior dependent on the initial condition in the range 
g^/n < T < Tg. But some discrepancies arise as r ap- 
proaches the lower bound of the temperature window. 
That is, as T approach Tg — Tsw in the left panel, (j) be- 
come smaller than predicted and tend to saturate. In 
the right panel, twin states were realized at the points 
A2 and B3, though they are in the lozenge phase in our 
theory. 

In Fig. 13, we show typical steady patterns of ip and 

(3i in the stripe states Ci, C2, and C3 taken at t/tQ = 
5 X 10'^. The patterns of 62 (not shown here) are almost 





FIG. 13: Stripe patterns of ^ (left) and ei (right) at the 
points Ci, C2, and C3, where t/tq — 0.6, —0.4, and —1.4, 
respectively, with a/ao — 1.4 and v/vo = 1 being common. 
In the left panels, ip is positive in the white regions, negative 
in the black regions, and nearly zero in the gray regions. Our 
theory indicates ■;/> = i/^s, —il^s, or in the three regions. In the 
right panels, ei is positive in the white regions and negative 
in the black regions. Theoretically, ei = e? + fig > or ei = 
Ci < 0. See the corresponding in the left panel of Fig. 12. 



indistinguishable from those of ip- The patterns in Fig. 
13 consist of long stripes for small (j>, but they contain 
even short ones and are considerably distorted for large (j). 
Thus, as T approaches Tg — Tsw, the domain morphology 
becomes more complex than in our theoretical picture. 
This explains why the data of (j) in the left panel of Fig. 
12 are below the theoretical straight line. 

In Fig. 14, the lozenge pattern at the point Ai is asym- 
metric with rectangular disordered regions, while that at 
the point B2 is more symmetric and some disordered re- 
gions are nearly square. The (f> for the pattern at B2 is 
0.51 and is larger than those at Ai and Bi in Fig. 12 
(being 0.42 and 0.43, respectively) by 20%. The shapes 
of the lozenges at B2 are more rounded than those at Ai 
and Bi with flattened edges. As a result, (j> can be larger 
at B2 than at Bi. 

We stress that the intermediate patterns strongly de- 



FIG. 14; Lozenge patterns of tp and ei at the points Ai and 
B2 in Fig. 8, where (r/ro, a/oo, d/uo) = (0.6,0.76,0.1) and 
(0.6, 1, 1), respectively. They are asymmetric ans symmetric, 
respectively, as in Fig. 5. At B2 the edges are flattened, 
resulting in an increase of the ordered regions. 

pend on how they were created. They are sensitive 
even to the initial noise. For example. Figs. 11 and 
14 show the lozenge patterns at the same point Ai at 
t/to = 5x10'^, where the large difference arises only from 
the initial noise. However, the pattern in Fig. 14 for Ai 
will change into a more symmetric one at longer times. 
Note that the lozenge angle 6 and the areal fraction (p are 
rather insensitive to the history. In fact, (f> = 0.41 and 
0.42 for the point Ai in Figs. 11 and 14, respectively. 



V. SUMMARY AND REMARKS 

Though in 2D, we have presented a theory of inter- 
mediate states at structural phase transition using a 
minimal model with a one-component order parameter 
-0 twofold coupled to the tetragonal strain properly 
(.9 7^ 0) ^iid to the dilational strain improperly (a ^ 0). 
The pinning of the domains is due to the nonlinearity 
(third and fourth order terms in the free energy) and 
there is no impurity. We summarize our main results. 

(i) In Eq.(2.27) or in Eq.(3.3) we have derived the elastic 
free energy contribution expressed in terms of ■0. It 
is very complicated due to the simultaneous presence of 
the proper coupling (oc g) and the improper coupling 
(oc a). 

(ii) For the stripe intermediate phase, the order parame- 
ter ip — ips within the ordered domains is determined by 
Eq.(3.15). The areal fraction 0, the transition reduced 
temperature Tg, and the width of the intermediate region 
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Tsw are expressed as in Eqs.(3.18)-(3.20). The free 
energy decrease is of the simple form Eq.(3.22). 

(iii) For the lozenge intermediate phase, the order 
parameter ip = ipi within the ordered domains is given 
by Eq.(3.38). The strains are approximately given by 
Eq.(3.35), where F in Eq.(3.30) is small within the 
ordered domains. The areal fraction 0, the transition 
reduced temperature r^, and the width of the interme- 
diate region T£^ are given in Eqs.(3.39)-(3.41). The free 
energy decrease is of the simple form Eq.(3.43). 

(iv) We have compared the free energies in Eqs.(3.22), 
(3.43), (3.48), and (3.52) for the intermediate and 
ordered phases in the plane of t and a with the other 
parameters held fixed. The transition reduced tem- 
perature, the phase diagrams, and the bird views of V' 
and 4> are in Figs. 7-9, respectively. The intermediate 
states can appear in the case a >aoi where ao is given 
in Eq.(2.30) and can be small near the tricritical point. 
Note that the ratio of the dilational strain and the 
tetragonal strains is of order fi/K for a ^ oq. 

(v) We have presented simulation results in 2D by 
integrating the dynamic equation (2.12). Figures 10 and 

11 illustrate the freezing processes of the domain growth 
resulting in mesoscopic intermediate states. Figure 

12 displays the areal fraction of the ordered regions 
obtained numerically. Figures 13 and 14 give typical 
patterns of the stripe and lozenge intermediate phases. 
These simulation results are consistent with our theory. 

Further we give some remarks. 

(1) In calculating the free energy F in the intermediate 
states, we have neglected the surface tension contribu- 
tion. This should be justified when the domain size much 
exceeds the surface thickness of order £ in Eq.(2.29). 

(2) The previous 2D simulation in the strain-only theory 
[U [TT] already produced the lozenge patterns in the 
presence of the third order term proportional to eie^, 
which plays the same role as the dilational improper 
coupling in this paper. 

(3) If a solid is compressed along one of the principal 
axes, the tetragonal extension along the compression 
axis becomes energetically unfavorable [21 [19], leading 
to a square-to-rectangle transition perpendicular to 
the axis. Then the transition could be described by 
a 2D theory with a one-component order parameter, 
although the lateral elastic deformations are much more 
complicated in real epitaxial films than treated in this 
paper [20] . 

(4) In its present form, our theory cannot explain the 3D 
experiments [31ll[51[Sl|71[Hl[n]- We should construct a 
3D theory of intermediate states. We will shortly report 



simulation results of intermediate states in 3D. 
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APPENDIX A: ORDER PARAMETER IN 
STRIPE INTERMEDIATE STATES 

We here solve the cubic equation (3.15). For u > and 
w = 0, we simply have 4's — 2ag/uL. Including the case 
of M < 0, we may generally solve it in the scaling form, 

V's = V't*(uM), (A.l) 

where we define two quantities, 

V-t = i3ag/2Lvy/^ (A.2) 
ut = 4z;^2/3 = (t;/3)i/3(4ag/L)2/3. (A.3) 

The scaling function ^(w) is determined by 

^{w)^ + w^{w) ^ 1 (A. 4) 



with $(?ii) > (since > 0). This cubic equation can 
be solved to give an explicit expression for <I>(w). For 
w > -3/4I/3 ^ -1.89, it reads 




For w < —3/4^/"^, it is expressed as 
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[1 -if^ 
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3/2y 
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3 V2 


w 





We find <i>(w) ^ l/w for w > 1 and <i>(w) ^ for 
w ^ —1, so that 

tps — 2ag/uL (w^wt), 

^ (3|«|/4t;)i/2 ((u«-Mt)- (A.7) 

The first line holds for u ^ Ut + /3l, while the second line 
for /3l > u + wt. 
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